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Abstract 



In this paper we develop a method for the estimation of Transport Induced by the 
Mean-Eddy interaction (TIME) in two-dimensional unsteady flows. The method is 
^^ built on the dynamical systems approach and can be viewed as a hybrid combination 

i-^ of Lagrangian and Eulerian methods. The (Eulerian) boundaries across which we 

I ^! consider (Lagrangian) transport are kinematically defined by appropriately chosen 

streamlines of the mean flow. By evaluating the impact of the mean-eddy interaction 
►j^ on transport, the TIME method can be used as a diagnostic tool for transport 

O processes that occur during a specified time interval along a specified boundary 

'^ segment. 

lO We introduce two types of TIME functions: one that quantifies the accumula- 

te tion of flow properties and another that measures the displacement of the transport 

geometry. The spatial geometry of transport is described by the so-called pseudo- 
lobes, and temporal evolution of transport by their dynamics. In the case where 
the TIME functions are evaluated along a separatrix, the pseudo-lobes have a re- 
lationship to the lobes of Lagrangian transport theory. In fact, one of the TIME 
functions is identical to the Melnikov function that is used to measure the distance, 
C^ at leading order in a small parameter, between the two invariant manifolds that de- 

flne the Lagrangian lobes. We contrast the similarities and differences between the 
TIME and Lagrangian lobe dynamics in detail. An application of the TIME method 
is carried out for inter-gyre transport in the wind-driven oceanic circulation model 
and a comparison with the Lagrangian transport theory is made. 
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1 Introduction and motivation 



Lagrangian transport methods are based on following the individual trajec- 
tories obtained by solving the original differential equation (ODE) for the 
particle location starting from a set of initial conditions Xq at time tg- 



— x = u(x,t) , (1) 

where u is the velocity field. The geometrical approach of dynamical sys- 
tems theory is particularly useful when the flow field has Lagrangian coherent 
structures that separate the flow into distinct regions. Then Lagrangian lobe 
dynamics describes the transport process between these regions using stable 
and unstable manifolds of hyperbolic trajectories as (moving) boundaries. The 
Lagrangian methods have been applied successfully to a number of unsteady 
geophysical flow problems; for a review, see [1,2]. If the flow is steady, i.e., 
u(x, t) = u(x), the invariant manifolds are stationary and no transport occurs 
between the regions. 

On the contrary, Eulerian-based methods are mainly concerned with the amount 
of transport across stationary (Eulerian) boundaries without computing indi- 
vidual trajectories. An advantage of Eulerian methods is that they tend to 
be much less elaborate than Lagrangian methods in terms of computational 
implementation. The choice for the Eulerian boundaries is generally flexible, 
unlike the Lagrangian methods. 

From the dynamical systems point of view, a parallel development of a method 
that computes transport across the Eulerian boundary has yet to take place. 
In this paper we begin the development of such a method. The method makes 
use of the interaction between the reference (mean) state and the unsteady 
variability (eddy) as the fundamental mechanism of transport. Hence we re- 
fer to it as the Transport Induced by the Mean-Eddy interaction (TIME). 
Using a streamline of the reference state as the boundary across which we 
consider transport, TIME can be thought as a hybrid of Lagrangian and Eu- 
lerian methods. Like the Eulerian method, the boundary is stationary. Like 
the Lagrangian method, the boundary is kinematically defined and there is no 
TIME in the steady flow without the unsteady eddy component in the veloc- 
ity. In certain situations we are able to describe the geometrical relationship 
of TIME along the Eulerian boundaries with Lagrangian lobe dynamics. 

We require no assumption of incompressibility in our theoretical framework. 
Therefore the ideas and techniques of the TIME method can be applied to two- 
dimensional compressible flow or three-dimensional volume-preserving flow 
which can be represented as special classes of two-dimensional flows, such as 



the shallow-water model. Remarks concerning incompressibility are provided 
throughout the paper as special cases. Extensions to three-dimensional flow 
are possible [3], but there is more complexity in the geometry of the transport, 
and this is will be the topic of a future publication. 

The outline of this paper is as follows. In Section 2, we provide a brief math- 
ematical background and introduce the notion of a kinematically-defined Eu- 
lerian boundary; readers who are familiar with elementary dynamical systems 
theory may omit this section without significant loss of continuity by refer- 
ring back to the notation and definitions as necessary. A brief glossary is also 
provided in Table 1. The TIME method is defined in Section 3, along with 
the two types of TIME functions. These functions, along with the notion of 
pseudo-lobes, are further explored in Section 4. An application of the TIME 
method is carried out in Section 5 for the inter-gyre transport in the double- 
gyre ocean circulation model and a comparison with Lagrangian transport 
theory is presented. Appendix A.l provides details of perturbation theory, 
and Appendix B compares the TIME method with the Lagrangian transport 
methods. [Tab.l] 

While in this paper we focus on introducing and developing the two TIME 
functions that estimate the amount and the geometry of transport, in the 
companion paper [4], we expand the TIME method further as a diagnostic 
tool for transport processes by analyzing in detail the infiuence of the mean- 
eddy interaction. 



2 Mathematical background 

In this section, we introduce the basic mathematical background necessary to 
develop the TIME method. The starting point is first expressing the velocity 
field (1) in the following form: 

u(x,t) = u(x) + u(x,t) (2) 

where {■} and {■}', respectively, correspond to the steady reference state and 
the unsteady fiuctuation around the reference state. The choice of reference 
state may not be unique. We choose the time-average (mean) of the full time- 
dependent field as the reference state in this study because the mean-eddy 
decomposition is natural when the fiow field is given by a data set; the TIME 
method itself does not require the reference state to be the mean. Many of our 
results will be perturbative in nature, with the (small) perturbation parameter 
being the amplitude of the fiuctuation that is implicitly included in u'(x, t) 
with respect to u(x). Appendix A gives results on the length of time intervals 



on which perturbed trajectories remain close to trajectories of the reference 
state. These results will provide the validity of the perturbative nature of 
our method since the TIME functions that we derive will be of the form of 
integrals along perturbed trajectories and approximations that, in principle, 
can be analytically computed are of the form of integrals along trajectories of 
the reference state. The regularity assumptions required on the velocity field 
are minimal. Essentially, we need existence and uniqueness of fluid particle 
trajectories, the ability to linearize about points in space, compute Taylor 
expansions through second order with respect to parameters, and for certain 
integrals of components of the velocity field along trajectories of the reference 
velocity field to exist. Assuming that the velocity field is twice continuously 
differentiable with respect to the spatial coordinates, time, and any parameters 
is adequate. No further assumptions on the nature of the time dependence (e.g. 
time periodicity, quasiperiodicity, etc.) are required. 

The use of perturbation theory in the development of the TIME method is 
made more transparent if we introduce an "order parameter", e, associated 
with the fluctuation term as follows: 



u'(x, t) = £ii(x, t) (3) 

where |u'(x'-^(s), t)| = (9(e|u(x'^(s))|). The introduction of e in this way makes 
perturbation arguments more transparent. However, the TIME functions can 
be equally as well expressed in terms of u'(x, t) or eu(x, t), but in either case 
the approximation is to leading order in the size of the fluctuation. 

2.1 Reference state and kinematically- defined Eulerian boundary 



We refer to a curve as Eulerian if it is stationary. The TIME method uses an 
Eulerian curve C = {x'"(s)} that is defined kinematically as a streamline of 
the reference flow. It can be given as a solution of 



^x^(s)=u(x^(s)) (4) 

with an initial condition xq = x'-^(so) at time to- For C to be a physically 
meaningful boundary, x'^(so) must be a regular point of n(xo), i.e., |n(xo)| 7^ 0. 

A trajectory with an initial condition x'^(so) at time to reaches x'^(so — to + t) 
at time t in the reference flow. This trajectory is uniquely identified by a scalar. 
So — to, because time shifts of a trajectory remain on the same trajectory in the 



reference flow. Throughout the paper, we interpret the flight-time coordinate 
variable s strictly as a spatial coordinate variable along C while t is a temporal 
variable. Accordingly (x'"(s),t) = (x'"(so — ^o + t),^) and 



{s,t) = {so-to + t,t) (5) 



can be viewed as different parametrizations of the same trajectory, which we 
call the reference trajectory. The Cartesian pair of coordinates {s,t) will prove 
to be particularly convenient for describing the TIME method. A glossary is 
provided in Table 1 for the principal definitions. 

A hyperbolic stagnation point is a singular point. In the reference flow, it is 
a special trajectory called the distinguished hyperbolic trajectory (DHT) [5] 
and we denote it by x'^''*, i.e., |il(x'^'^*)| = 0. Although x^'^* itself cannot be a 
physically meaningful Eulerian boundary C, the unstable and stable invariant 
manifolds that have a DHT at the starting and end point, respectively, are 
special types of C: 



W^ = {x^(s)| lim x^(s) = x'l^n; (6a) 

s— >— oo 

W^ = {x^(s)| Jinn x^(s) = x^^^*}, (6b) 

where W denotes the reference manifold with the superscripts {■}'^ and {■}^ 
for unstable and stable invariant manifolds, respectively. In a case where W 

g 

and W coincide, the reference streamline is called the separatrix or hetero- 
clinic connection of the upstream DHT x^l^ and the downstream DHT xj^*: 

W"" = {x^(s)| Jm^x^(s) = x'i'j^,]!™^ x^(s) = x^} , (7) 

where the superscript {■}^ stands for heteroclinic connection and the sub- 
scripts {-jipoo represent the direction of s towards the corresponding DHT. In 
addition, if xl'j^ and x^* coincide, then W is called the homoclinic connec- 
tion. An invariant manifold is special case of kinematically-defined C because 
s has a semi-infinite or bi-infinite range as in (6) and (7). We emphasize that 
the terms finite and infinite refer to the range of s on C, rather than the 
physical length of C. 

For the description of the transport geometry near C, it is often convenient 
to use an orthogonal arc-length coordinate system, (/,r). Along C, the arc- 
length / = l'"{s) and the flight-time s are related by the local velocity, i.e.. 



^/*"(s) = |u(x'^(s))|. Normal to C, r is defined to be the signed distance of a 
neigliboring point x to C; r > 0, r = and r < correspond to tlie left, on, 
and the right of C with respect to the forward direction of x*" (s) along C . A 
pair of orthogonal unit vectors in the tangent and normal directions to C are 
given by 






where u = {ui,U2)'^ . The transformation between the Cartesian and arc-length 
coordinates is area-preserving. 

2.2 Unsteady flow and perturbation theory 



As we have noted, the mathematical formulation of the TIME method is based 
on perturbation theory for a velocity field given by (2) and (3). The necessary 
background and results are given in Appendix A. Trajectories of the unsteady 
flow passing through xq on C at tg are of the following form: 

x(t; x^(so), to; e) = x^(so - to + t) + £x(t; sq, to) + 0{e^), (9) 



where £:x(t; So,to) is the leading-order displacement vector with x(to; Soj^o) = 
0. Computing the Taylor expansion of (2) and the time derivative of (9) with 
respect to e gives to the following linear ordinary differential equation for 
x(t;so,to): 

-x(t; So, to) = ^xu (x^(so - to + t)) X + u(x^(so - to + t), t). (10) 

Given C = {x'^(s)}, x(t; So,to) can be obtained by solving this linear system 
where the nonlinear evolution of x'^(so — to + t) provides us with the time- 
dependent coefficients and the inhomogeneous term. In Appendix A we show 
that perturbation theory can provide valid approximations in situations where 
C is defined over finite, semi-infinite or bi-infinite time intervals. 

It is worth noting here that for many of the most fruitful perturbation theories 
used in dynamical systems type analyses rarely are precise bounds available for 
the size of the perturbation for which the method is applicable. Nevertheless, 
this has not limited the insights they have provided in a variety of applications. 
For example, the typical statements of Melnikov's method [6,7,8] indicate only 



that it is valid for e sufficiently small. Another example is the well-known 
Kolmogorov-Arnold-Moser (KAM) theorem [9], which has been proven useful 
in many applications despite the fact that the bounds are generally too strict to 
be practically applicable. The situation with the KAM theorem is even worse 
since rarely are the hypotheses of the theorem even verified in applications 
since they require the velocity to be expressed in action-angle variables, which 
can rarely be achieved. This limitation also prevents one from obtaining any 
type of bound on the perturbation for which the theorem is valid. 



3 Transport functions for TIME 



Having kinematically defined the Eulerian boundary C by the reference state, 
we now turn our attention to transport across C . There are two aspects: one is 
concerned with the amount of flow property and the other is concerned with 
geometry of transport. Examples of flow properties are mass, temperature, 
humidity in the atmosphere, salinity in the oceans, and such. The TIME func- 
tions are developed for these two aspects, ffist for a finite time interval along 
any C (Section 3.1) and then for an infinite time interval along an infinite C 
(Section 3.2). 



3. 1 Derivation of the finite-time TIME functions 



3.1.1 Accumulation of a flow property 

We assume that the time-dependent fluctuation in the flow property distribu- 
tion, denoted by (^(x, t), is also small 



Q(x,t) = Q(x) + Q'(x,t). (11) 



Like u(x, t) in (2) and (3), we introduce Q'(x, t) = e(5(x, t) with |(5'(x'^(s), t)| = 
0{e\Q{^^{s))\). To illustare the basic idea for estimating the amount of prop- 
erty transport, we consider the imaginary fluid column J-" in the flow (Fig- 
ure Ic). By accumulating the flux at the moving intersection of J-" with C, we 
obtain the net amount of accumulation. [Fig.l] 

Up to leading order, the intersection (s, t) of J-" with C at time t is approx- 
imated by the reference trajectory (sq — t^ + t, i) using perturbation theory 
(Appendix A.l). At (s, t), the instantaneous flux of (5(x'"(s), t) carried by the 
local velocity u(x'-^(s),t) across C per unit length is 



ef(/^(s)) A [ g(x^(s)) 5u(x^(s),t) +e2g(x^(s),t) u(x^(s),t) ] , (12) 



where $,±{l'"{s)) and ^n (/'"(s)) are defined in (8). The positive value means the 
fiux from the right to the left across C with respect to the forward direction 
of s. This formula (12) says that the instantaneous fiux of Q across C exists if 
u(x, t) has a component normal to C, and that the time-dependent fiuctuation 
Q\'x.,t) = eQ{'x,t) contributes to the transport at the higher order. At the 
leading order, the instantaneous fiux of Q penetrating across C at (s, t) per 
unit fiight time is g(x'^(s)) n(x'^(s)) A £u(x'^(s), t), i.e., 



Q(x^(s)) u(x^(s)) A u'(x^(s), t) = g^(s) /i^(s, t) (13a) 
using (3), where 

g^(s)=Q(x^(s)) (13b) 

/i^(s,t)=u(x^(s)) Au'(x^(s),t). (13c) 



We refer to jjp{s,t) as the instantaneous flux function, induced by the un- 
steadiness (eddy) of the velocity through the interaction with the reference 
(mean) fiow. This is the origin of the transport induced by the mean-eddy 
interaction (TIME) across C. The sign of fi'"{s,t) indicates the direction of 
the instantaneous fiux. 

The accumulation over the interval [to, ti] is thus approximated by //^^ q'"{so — 
to + t) /i*" (sq ^ ^0 + t, T)dT up to leading order. Because this amount is the 
same for any (s, t) along the reference trajectory with s — t = SQ — to,we obtain 
a general form of the accumulation 



m'^'{s,t;to:ti)= q^{s - t + r) fx'^\s - t + T,T)dT , (14) 

Jto 

where the first pair (s, t) in the arguments of the left-hand side represents 
the combination of spatial coordinate and time at which the net accumulation 
of Q is evaluated, while the next pair (to : ti) concerns the time interval on 
which the transport takes place. Here t can be either before or after to or ti. 
We refer to m'"(s,t;to : ti) as the accumulation function. Characteristics of 
m*" (s, t; to : ti) will be discussed further in Section 4. 
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3.1.2 Displacement distance and area 

For the geometry, we consider the displacement distance of the particle starting 
from {so,tQ) on C. In the unsteady flow at time t, the displacement of the 
particle from x'-^(so — ^o + ^) is x'(t; x'^(so), ^o) up to leading order by (10). For 
particle transport and its geometry, we choose to use arc-length coordinates 
(/, r) in the description of the displacement functions because the displacement 
distance has the physical dimension of length. Using (9) along a reference 
trajectory (s,t) = (sq — ^o + ^) with the initial condition (/qj^o) = (^'"(■5o))^o)) 
the leading order term for the displacement distance due to particle transport 
at t G [toj^i] is given by 



r'{t; /o, to) = If (/^(so -to + r)) A x'(t; sq, to) ■ (15a) 

Because x'(t; So, to) = £:x(t; sq, to) can be obtained by solving (10), then so can 
r'(r;/o,to) by the direct substitution. However, a simpler formula is available 
by considering 



a(t; lo, to) = r'(t; Iq, to) |u(x'^(s)) | = u(x'^(so - to + t)) A x (t; sq, to) , 

(15b) 

which corresponds to the displacement area per unit s along C as shown in 
Appendix A. Using a' (to; /o^to) = for the initial condition, construction, the 
solution for the displacement area at ti is given by 



a'(ti;/o,to)= / e^'{s-t + ti: s-t + r) fi'^{s-t + T,T)dT , (16) 



where 

e^{s : So) =exp{ /' trace{i:'xn(x^(^))}rf^} (17) 

JSQ 

reflects the compressibility of the reference flow; for an incompressible flow, 
e^(s : So) = 1- 

Like the accumulation, the displacement area associated with transport over 
[to,ti] can be evaluated at (/'^(s),t) where t can be before, in, or after the 
time interval. Conceptually, this is to let a'(ti;/o,to) obtained by (16) evolve 
under the reference flow over an additional time interval [ti,t] to take the 
incompressibility into account. As shown in Appendix A, the flnal form of the 
displacement functions is given by 
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a'^{s, t; to : ti) = e'^(s : Sq - to + ti)a'{ti; Iq, to) 

= / e'^(s : s-t + r) /i^(s - t + t, T)dT, 

Jto 

(18a) 



C/iC/ \ ^ ^ ^ \ _'^ {s,t;to : ti] 

|u(x^(.))| 



r^(/^(s), t; to : ti) = .r.L-c " \i ' (l^^) 



As in the case of m'"(s,t;to : ti), the first pair (/'^(s),t) in the argument 
represents the spatial coordinate and time at which the function is evaluated, 
and the next pair correspond to the time interval when transport takes place. 

Accordingly over [ti,to], the displacement is determined by two contributions: 
one is from the unsteadiness of the flow measured along C through the in- 
stantaneous flux fi'"{s — t + t,t), and the other from the compressibility of 
the reference flow through e'"(s : s — t + r), which may result in compression 
or expansion of the area. Sign of a'-^{s, t; to : ti) indicates the directionality of 
transport across C. 



3.2 Extension over the infinite TIME functions 



We refer to the accumulation function (14) and displacement distance and 
displacement area functions (18) as the (finite) TIME functions because they 
are defined over a finite time interval and hence the finite range of s along C. 

If g 

These TIME functions can be extended over the semi-infinite W and W and 
bi-infinite W , because the exponential decay of the velocity towards DHTs 
at the starting or (and) end point (s) guarantees the convergence conditions 
required for the validity of perturbation theory (Appendix A). The extension 

TJ 

of the displacement functions over W can be particularly useful since it pro- 
vides a direct link to the Melnikov function which measures the leading order 
distance between the time-dependent unstable and stable invariant manifolds 
[6,7,8]. The Melnikov function has been used to study Lagrangian transport, 
mostly in incompressible flows [10]; also see [2]. 



The TIME functions for these special C are as follows. The transport that 
has happened in the past across W can be obtain 
functions over a semi-inflnite time interval (— oo,tn 



has happened in the past across W can be obtained by extending the TIME 
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m^{s,t;to)= q'^{s-t + T) n^\s-t + T,T)dT; (19a) 



a^(s. 



t; to) = / " e^{s : s-t + r) ff {s - t + t, T)dT; (19b) 

for some tg- Similarly, the transport that will happen in the future across the 
stable manifold W can be obtained by extending the TIME functions over 
the semi-infinite interval [to,c>o): 



m (s. 



t;to)= -f {s-t + r) ^^{s-t + T,T)dT; (20a) 



a''{s,t;to)= e^{s : s -t + r) iJ,^{s -t + T,T)dT; (20b) 

Jto 

Finally, the entire transport across the separatrix W can be obtained by 
extending the TIME functions over the bi-infinite interval (— oo, oo): 

m"{s, t) = m^{s, t; to) + m^{s, t; to) 

/oo 
q^{s-t + T) iJ,'^{s-t + T,T)dT; (21a) 

-oo 

a"{s, t) = a^{s, t; to) + a'^(s, t; to) 

/oo 
e'^\s : s-t + r) iJ,'^{s-t + T,T)dT; (21b) 

-oo 

r^{l^{s),t) ^ T^l^{s),t; to) + r^(P(s), t; to) = SS^- (21c) 



The displacement area function a^{l^{s),t) for W is the same as the so- 
called Melnikov function. 



4 Characteristics of TIME 



4-1 Characteristics along an individual trajectory 



The accumulation is obtained by following the individual reference trajectories 
(Section 3). This leads to the concepts of invriance and piece-wise indepen- 
dence. 
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Invariance of the accumulation function. Given a fixed time interval [to,ti], the 
accumulation function is invariant: 



m^ (s, t; to : ti) = m^ (s - t + r, r; to : h) (22) 

for any r. This invariance implies that each trajectory has perfect memory for 
the amount of transport. Invariance for the displacement area function: 



a^{s, t;to : h) =e^ {s : s - t + t) a'^ {s - t + r, r; to : ti) (23) 

is subject to the compressibility factor e^'{s : s — t + r) of the reference flow. 

Piece-wise independence. It is clear in the deflnition (14) that the time interval 
[to, ti] can be broken up into an arbitrary number (say K) of pieces: 

K 



m [s 



t; to : ti) = ^ m^{s, t; to.k-i : to.k) , (24) 



k=l 



where U^;^[to.fe-i,to.fc] = [to,ti] with to.o = to and to._ft: = ti. This is a temporal 
piece-wise independence. Using (5) along {sq — to + t, r), spatial piece-wise 
independence follows naturally by breaking the spatial segment [sa, Sf,] into 
K pieces by uf^Jsa.fc-i, Sa.fc] = [sa,Sb] with Sa.o = Sa and Sa.K = Sb and 
transforming them into K temporal pieces with [to.fc-i, to.fe] = [to — Sa.o + 
Sa.k^i, to - Sa.o + ^a.fe] ovcr [to, ti] = [to, to - Sa + Sfo]. Bccause the independence 
is a characteristics deflned for a fixed {s,t), both temporal independence and 
spatial independence hold for the displacement area function a'"(s,t;to : ti). 



4-2 Coherency of transport 



Geometry of particle displacement leads to the concepts associated with the 
coherency of transport. Because the geometrical characteristics discussed here 
hold for any C over any time interval, we will drop [to : ti] from each notation 
for simplicity after the first appearance as indicated in [= {■}]; for example, 
r^{l'^\s),t)[= r'^(/'^(s),t;to : ti)]. When C is taken as W , there is a geo- 
metrical relation to Lagrangian transport, which we will treat separately in 
Appendix B. 



For an illustration of transport geometry and coherency, let us consider an 
imaginary material curve TZit) placed initially on C at time to, i.e., TZito) = 
{{l,r)\r = 0} (Figure Id). In the reference fiow, TZ{t) advects along C without 
any displacement. In the unsteady flow, velocity normal to C may let 7^(t) 



14 



depart from C . For the transport geometry associated with the TIME method, 
we define 



i?^(t) = {(/, T)\l = f{s),T = r^(/^(s), t)} , (25) 

which is the leading order approximation to TZ{t). 

Pseudo-primary intersection point (pseudo-PIP) sequence of R'"{t) and C. In 
the unsteady fiow, R^{t) may intersect with C to form a chain of lobe-like 
structures as shown schematically in Figure Id. We denote such a ordered 
sequence of such zeros {s'^{t)[= s'^it] to : ti)]} with 

sf (t) = {s I a^(s,t) = 0, sf (t) < s$+i(t)} , (26) 

where we use the fact that the zeros of a'"{s^ t)[= a^{s, t; tg : ti)] and those of 
r^{l,t) are identical. Unless the zero is non-degenerate, {s^(t)} is generally 
identical to the intersection sequence of TZit) up to the leading order; see [11] 
and also Appendix A.l. We call {s^(t)} the pseudo -primary intersection point 
(pseudo-PIP) sequence in contrast to the Lagrangian lobe dynamics for the 
heteroclinic connection. The term 'primary' is used here to emphasize the 
analogy of the PIP of Lagrangian lobe counterpart. A pseudo-PIP sequence 
can be transformed into arc-length coordinate /^(t)[= /^(t;to : ^i)] so that 
/f (t) ^ /^(sf (t)). [Fig.2] 

Invariance of pseudo-PIP sequence. Because the displacement area a'"{s,t) is 
invariant subject to the compressibility effect from (23) and the compressibility 
effect will not change {s'^ if)} at given t, {s^{t)} is invariant. Each s*^ {t) 
coincide with a reference trajectory, i.e., 

(5f(r),r) = (sf(t)-t + r,r). (27) 



Pseudo-lobe seqence. We denote the lobe-like structure defined by the segments 



of R'"{t) and C between a pair of two adjacent pseudo-PIPs by L^j-^^(t)[E 
L9j+i(t;to : ti)] and call it the pseudo-lobe: 



LlAt) = {{l,r) I r{r-r^{l,t)) < 0, / G [/f (t), /f+,(t)]} . (28) 

Its sequence makes the chain-like structure, which we call the pseudo-lobe 
sqeuence, {L^j^iit)}. Using a'~"\s,t), the pseudo-lobe is the area surrounded 
by a*"(s,t) and C between two adjacent pseudo-PIPs. It is worth emphasizing 
that temporal and spatial piece-wise independence (Section 4.1) holds for the 
displacement function of each pseudo-lobe individually. 
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Signed area. The size of each pseudo-lobe measures the amount of the locally 
coherent transport. Its signed area is given by: 



/•«'?,, (t) /-s*?, , (t) 

If the flow is incompressible, the area of each lobe is invariant, i.e., A(L^-^]^(r)) = 



^(^^•+iW). 



Directionality of pseudo-lobes. Each pseudo-lobe represents the amount of fluid 
particles that go across C over [tc^i]- Depending on whether R^{t) lies to 
the left or right of C, L'^-j^^it) can be of two types, denoted respectively by 
j^RC Qj, j^CR corresponding to r^{l,t) > or r'-'\l,t) < for transport across 
C from right to left or left to right. Provided all intersections are transverse, 
{L^j_,_]^(t)} alternated the types between L^*" and L*"^ along C, resulting in 
a chain structure at any given time. 

The geometry associated with transport of Q can be established concep- 
tually by replacing the displacement area function a^{s,t) for R'"{t) with 
the accumulation function m'"(s,t)[= mP{s,t;tQ : ti)] in arc-length coordi- 
nates. In the case for mass with g'"(s) = 1 where the flow is incompressible, 
this operation results exactly in the displacement functions of fluid particles, 
rrf{s,t) = a'^{s,t). 



Application to a numerical simulation of the wind-driven double- 
gyre ocean circulation 



In this study we apply the TIME method to the inter-gyre transport in the 
mid-latitude, wind-driven ocean circulation. The data set is obtained by a 
numerical simulation of a quasi-geostrophic (QG) 3-layer model in a rect- 
angular basin geometry with free slip boundary conditions [12]. Due to the 
latitudinal antisymmetric wind-stress curl applied at the ocean surface, the 
basic circulation pattern in the top layer is a double-gyre structure separated 
by an eastward jet shooting off from the confluence point of the southward 
and northward western boundary currents (Figure 3). Driven by this strong 
jet, the subpolar gyre circulates counterclockwise in the north and the sub- 
tropical gyre circulates clockwise in the south. Depending on the value of the 
parameters such as the viscosity and the wind stress curl, the ocean circulation 
exhibits a rich time-dependent dynamics [13,14]. [Fig-3] 

At a wind stress curl of 0.165 dyn/cm^, the ocean dynamics is nearly periodic 
with dominant spectral peak at period T ^ ISldays in a 1000 km x 2000 
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km rectangular domain. We choose this flow since the physical interpretation 
for a flow fleld close to periodic is more simple and therefore it allows us to 
focus more on demonstrating our method. It is worth noting again that the 
method itself does not require time periodity of the flow. The velocity data 
set used in this study has a spatial resolution of 12. 5kmx 12.5km and is saved 
daily after a 30,000-day spin-up from rest. Figure 3a shows the streamfunction 
■?/'(x, t) of the top layer at t = ^944 when '?/'(x, t) is close to the reference state 
ipi'x). From here on, the subscript in t* denotes days after the completion of 
the spin-up. The most signiflcant region of unsteadiness of the flow lies in the 
upstream region of the eastward jet near the western boundary, while small- 
amplitude Rossby waves propagate westward in the entire ocean basin. The 
Lagrangian transport processes between the two gyres are governed by the 
lobe dynamics associated with the unstable invariant manifold W^(t) of the 
upstream DHT on the western boundary and the stable invariant manifold 
yV^{t) of the downstream DHT on the eastern boundary; see Appendices A.l 
and B for the deflnitions and more details of VV'^(t) and W'^(t). The inter-gyre 
transport in the top layer was carefully studied by [15] using Lagrangian lobe 
dynamics methods. 

For the comparison with the Lagrangian method, we choose the Eulerian 
boundary of the TIME method as the bi-inflnite, reference heteroclinic con- 
nection W = {x^(s)} which spans over s G (—00, 00) with x'^(— 00) = x'l'j^ 
on the western boundary and x^(oo) = x^^ on the eastern boundary (solid 
line in Figure 3c). It is worth noting again that application of the TIME 
method is not limited to flows that possess a heteroclinic connection. For the 
computation of the TIME functions, we choose the location of s = so that 
x^(0) is very close to x"!^ and |ii(x-^(0))| is exponentially close to zero. Along 

W , the speed |il(x-^(s))| of the mean jet signiflcantly increases starting near 
x^(sj). Geographically x^(sj) is separated from x-'^(O) only by 600m. As s in- 

TT 

creases towards the downstream direction, W makes a sharp turn in around 
x^(sn) in the north, followed by the second sharp turn around x-^(ss) in the 
south. Measured in the flight time, sj, sn, and ss are llOdays, 129days, and 

TT 

174.5days, respectively. In the further downstream direction, W extends to 
the east and exhibits little meandering. After about x-^(s25o) with the sub- 
script in s* for the flight-time coordinate in days from s = 0, |u(x-'^(s))| be- 

TT 

comes extremely small. In the region near W , the order of the unsteadiness 
relative to the reference state is small (~ 0.1), supporting the applicability of 
the TIME method (see Appendix A. 2). 

Figure 4 shows the instantaneous flux fJ^^{s, t) on W^ as a Hovmoller diagram 

[16,17] in the (s, t) space for [0, 500] x [750, 1250]. The signals of fi^i-s, t) are [Fig.4] 

periodic in t with the period T because of the time-periodic ocean dynamics. 

For small s, fi^{s,t) is always near zero because u(x^(s)) is exponentially 

small near xi'j^. Also for large s, n^{s,t) is always almost zero because both 
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u(x^(s)) and u'(x-^(s), t) are small there. Most signals of /x^(s, t) are confined 

in between. Figure 5 shows the four different phases of n^{s,t) during one 

period of the ocean oscillation at t = t^^^, tgg2(= tg44 + T/4), tio2o(~ ^^944 + 

T/2), and tio58(= ^944 + 3T/4). Downstream propagation of the large signals is [Fig. 5] 

seen over [sj, sn] where typically three extremes (one local minimum and two 

local maxima, or vice versa, depending on the phase of the ocean dynamics 

during T) propagate downstream. Over the subsequent segment, upstream 

propagation of the small signals is observed. These signals are a consequence of 

the double-gyre ocean dynamics (variability) through the interaction with the 

reference (mean) flow. In the companion paper [4], we examine the details of 

the variability in the double-gyre ocean dynamics and its relation to transport 

processes, while the focus of this paper is on the comparison of the results 

of the TIME method with those obtained by the Lagrangian lobe dynamics 

method. For the computation of the TIME functions, it suffices to use fi^{s, t) 

over a short finite segment rather than the entire (—00, 00) along W because 

fi^{s,t) decays to zero for both small and large s. Our results presented here 

use n^{s,t) over [sq, 5*000]- ^ straight diagonal line in Figure 4 is an example 

of the reference trajectory along which the TIME functions are integrated. 

Figure 6 shows the infinite displacement area functions. By choosing to = t [Fig. 6] 
in (19)-(21), a'^{s,t;t) is the transport that has already happened in the past 
of the present time t, a^{s,t; t) is the transport to happen in the future of t, 
and a^{s,t) is the net inter-gyre transport over all time. The two properties 
of the TIME functions (Section 4) are observed in this figure as follows. The 
temporal piece-wise independence property (24) is observed by a^{s,t) = 
a^{s,t;t) + a^{s,t;t) for any (s,t). The invariance property (23) is observed 
by a^{s, t) = a^ [s + 5,t + 5) for any 5. 

Due to the active fi^{s,t) over the segment [sj,ss], the displacement area 
functions have the relation a^{s,t) ~ a^{s,t;t) for s < sj. This is because 
no transport has happened there yet, leading to a'^{s,t;t) ^ 0. Similarly, 
a^{s,t) ^ a'^{s,t;t) holds for s > sg because all transport has occurred al- 
ready, leading to a''^{s,t;t) ~ 0. Because the bi-infinite displacement area 
function is the same as the Melnikov function, an implication of these results 
to the Lagrangian transport is that the development of the Lagrangian lobes 
are mainly governed by the fiow dynamics between x'^(sj) and x^(ss). The 
Lagrangian lobes simply advect in the downstream direction of x^(ss) without 
any further inter-gyre transport activities. 

Figure 7 shows the corresponding displacement distance functions, r^{l,t;t), 
r^{l,t;t), and r^{l,t), along the arc-length coordinate / = l^{s) measured in 
km. The displacement distance functions are inversely proportional to the local 
reference velocity |u(x-'^(s))|. Therefore, the pseudo-lobes defined by the two 
consecutive zeros of the displacement distance functions (see Appendix B) are 
stretched and stay near zero around s^^g where |u(x'^(s))| is maximum along 
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W . Accordingly it is near s = s^^g that the inverse pseudo-PIP s^{t) associ- 
ated with r'-' {l^ (s) , t) and r^{l^{s),t) exists for the pseudo-turnstile mecha- 
nism (see Appendix B). This also implies that the turnstile mechanism of the 
Lagrangian lobe dynamics for the inter-gyre transport occurs near x'^(s^;^g). In 
the downstream for large s, all displacement distance functions show vertically 
elongated pseudo-lobes due to very small |il(x^(s))|. The same phenomenon 
occurs in the upstream for small s near the upstream DHT. This is related to 
the elongated structures of the Lagrangian lobes near the DHTs (FigureSb). [Fig. 7] 

The pseudo-lobes of the bi-infinite displacement function can be used to es- 
timate the amount of the Lagrangian inter-gyre transport carried by the La- 
grangian lobes (Appendix B). Table 2 summarizes a quantitative comparison. 
The area A{Cj.j_^_l(t)) of the Lagrangian lobe CJJ_^_-^^(t) is computed by the La- 
grangian method [15], while the area A{L^,j_^_^(t))) of the pseudo-lobe L^,jj^i{t) 
is computed by (29). Agreement between the two methods is quite good. [Tab. 2] 

Clearly the Lagrangian lobe dynamics method provides precise geometry as- 
sociated with particle transport that the TIME method cannot (compare Fig- 
ure 3b with Figure 7). The computation of W^(t) and W^(t) using a velocity 
data set given on a grid can be computationally intensive. As described in 
[18], it requires sophisticated spatial and temporal interpolation schemes. In 
Figure 3b, VV'^(t) and VV'^(t) have about 5000 and 2000 particles, respectively. 
With this relatively large number of particles, it is also computationally chal- 
lenging to preserve the area of each Lagrangian lobe precisely, due to the 
geometrical complexity of W'^(t) and VV'^(t) as well as the limited accuracy of 
the numerical schemes. In contrast, the TIME method is extremely efficient 
in computing the transport because it requires a very small number of the 

TT 

simple operations. In this example, we used only 1001 data points along W 
to compute a^{s,t). Once a^{s,t) is obtained, it automatically provides the 
value a^{s + t — t,T) for any (s + r — t, r). 



6 Summary and concluding remarks 



We have developed a mathematical framework for the estimation of Transport 
Induced by the Mean-Eddy interaction (TIME) for flow properties and fluid 
particles with emphasis on two-dimensional unsteady geophysical flows, with- 
out the assumption of incompressibility. The TIME method estimates the 
amount of Lagrangian transport across the kinematically-defined Eulerian 
boundary by the appropriately chosen streamline of the reference flow (Sec- 
tion 2). The TIME method is a hybrid combination of Lagrangian and Eulerian 
methods and is based on the dynamical systems approach. It enables on to 
analyze unique features of transport that neither the Lagrangian nor Eulerian 
methods can provide. 
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By considering two different aspects of transport, we obtain the accumulation 
function for flow properties, as well as the displacement distance and area 
functions for fluid particles (Section 3). The dynamical systems approach leads 
to the useful characteristics such as invariance, independence and coherency of 
the geometry (Section 4). In the companion paper [4], we develop a framework 
for the analysis of the transport process in which these characteristics play a 
key role. 

The notion of pseudo-lobes is developed to describe the geometry associated 
with TIME. When a heteroclinic is used as the Eulerian boundary and the 
time interval for transport to take place is set over a bi-infinite time interval, 
the pseudo-lobes are geometrically closely related to the Lagrangian lobes of 
the associated invariant manifolds. The novel turnstile mechanisms for La- 
grangian transport can be carried over in the TIME method by taking the 
mirror image of the pseudo-lobes in the upstream region of the heteroclinic 
connection (Appendix B). An application to an oceanic problem and a com- 
parison with the Lagrangian lobe dynamics studied by [15] (Section 5). 

The TIME method is designed to augment and supplement Lagrangian and 
Eulerian transport methods by providing the unique capability to analyze 
the underlying transport processes, as it will be shown in [4]. The method 
can be applied to more genercal cases than the application presented in this 
paper; the method itself does not require the time periodicity of the flow filed 
or a heteroclinic connection in the reference state. Various applications and 
extensions of the Eulerian transport theory, including three-dimensionality [3] , 
brings a new point of view and direction to transport studies in geophysical 
flows. 
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A Mathematical Background on Perturbation of Trajectories 



In this appendix we given the necessary mathematical background that the 
TIME functions are valid approximations to the quantities that they mea- 
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sure to leading order in the size of the fluctuation around the reference state. 
Essentially, the result that we need is that trajectories of the full, time de- 
pendent velocity field are "close", in the sense of the size of the fluctuation, 
to trajectories of the reference state velocity field on the time intervals of in- 
terest. The TIME functions are integrals of functions involving the pieces of 
the decomposed velocity field, flow properties, and geometrical features of the 
chosen Eulerian curve C 

Two distinct situations need to be considered. One is where the velocity fleld 
is deflned by an analytical formula, discussed in Section A.l. The other is 
where the velocity fleld is deflned as a data set, discussed in Section A. 2. 



A.l A Velocity Field Defined by an Analytical Formula 



In this appendix we collect together the results that we use on the approxima- 
tion of trajectories of the reference flow by trajectories of the flow consisting 
of the reference flow and the fluctuations about the reference flow over appro- 
priate time intervals. These results are stated here for completeness and they 
can be found in a number of references. See, e.g., [19,20,7,8,21,22,6]. These re- 
sults deal with the case of time-periodic fluctuations. However, the arguments 
and proofs for aperiodically time-dependent fluctuations are the same and are 
discussed in [23,2]. 

First, we recall notation established in Section 2. 

The velocity fleld expressed analytically as the sum of a steady reference state 
and an unsteady fluctuation is given by: 



— X = u(x, t) = u(x) + £u(x, t), (30) 



(note that the following perturbation results are valid for x either two or three 
dimensional). We denote a trajectory of (30) by: 



x(t;x^(so),to). (31) 

The reference velocity fleld is given by: 

^x^(s)=u(x^(s)), (32) 
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We denote a trajectory of the reference velocity field by: 



x^(so-to + t). (33) 

Basic results in the theory of ordinary differential equations say that if (30) is 
a C" (i.e. r times continuously differentiable) function of t, x, and e, then (31) 
is a C^ function of t, s, and e. Hence, we can Taylor expand in any of those 
variables. 

In particular, (31) can be substituted into (2) and differentiated with respect 
to e. In this way we obtain ordinary differential equations for the coefficients 
(which are functions of time) for the different powers of e. Following this 
procedure, we obtain: 



x(t; x^(so), to) = x^(so -to + t)+ £x(t; Sq, to) + ^(e'), (34) 

where x(t; So,to) satisfies the first variational equation: 



-X = D,u(x^(so - to + t))x + u(x^(so -to + t),t), (35) 

which is a linear, inhomogeneous, differential equation. It is clear that with 
more work one could derive differential equations whose solutions are the 
coefficients of the higher order terms in the e expansion of of the trajectories. 
However, the first order term in e will be sufficient for our purposes. 

In order for our results to be valid we will need estimates on the time interval 
for which (31) and (33) are 0{e) close (as measured in an appropriate norm, 
say the maximum of the Euclidean distance between the two trajectories over 
the time interval of interest). 

We state our first result on closeness over finite time intervals. 

Proposition A.l (Finite Time Approximation.) Suppose \:si{tQ;x'-''{sQ),tQ 
x^(so)| = 0{e). Then |x(t;x^(so), to) - x^(so - to + t)| = 0{e) for \t -to\ = 
Oil). 

This result is elementary and well-known, and is a consequence of a simple 
application of Gronwall's inequality that can be found in the references given 
above, as well as many texts dealing with perturbation results. 

In Section 2 we discussed the situation where the (steady) reference fiow con- 
tained a hyperbolic stagnation point, denoted x'^'^*, with the hyperbolic stag- 
nation point having stable and unstable manifolds, W (x'^'^*) and W fx'^ 
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respectively. We next state a result that describes how this hyperbolic struc- 
ture persists under perturbation of the reference flow by the fluctuation. 

Let Beg denote the ball of radius eq centered at x'^'^*. Let W^^^ (x.'^^'^j denote 
the component of the intersection of W f x'^'^*) with i^^,, that contains x'^'^*. 
Similarly for W^^^ (x'^'^*]. Then we have the following well-known result. 



Proposition A. 2 (Persistence of Hyperbolic Structures) There exists Eq 
sufficiently small such that for all < e < Eq (30) has a hyperbolic (time de- 
pendent) trajectory x^'^*(t) = x*^"^* + 0{e). Moreover, x^'^'(t) has local stable 
and unstable manifolds, denoted W\^^ Vx-^^{t)\ and W^^^ (x^'^*(t) j, that are C" 

E close to W^ioc (x'^'^M and W^^^ fx'^'^^M, respectively. 

Note that it is an "infinite time" result. 

Using an argument that combines Proposition A.l and Proposition A. 2, the 
following result can be proven. 

Proposition A. 3 (Approximation on Semi-Infinite Time Intervals) Suppose 
we choose initial conditions of trajectories in the stable and unstable manifolds 
of the hyperbolic trajectory at t = t^: 



x^(t;x^(so),to)GW-''(xfHto)), (36) 

x^(t;x^(so),to)GlF''(xfHto)) (37) 



then 



x^(t; x^(so), to) = x^(so - to + t) + £x^(t; So, to) + 0{e^) 
uniformly in t for t E [to, oo) 

x^(t; x^(so), to) = x^(so - to + t) + £x^(t; So, to) + 0{e^) 
uniformly in t for t G (— oo,to] 



(3J 



where x'^(t; Sq, to) and x'^(t; So, to) satisfy 
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-x^ = Dxu(x^(so-to + t))x^ + u(x^(so-to + t),t) for te[to,oo), 

-x^ = Dxu(x^(so-to + t))x^ + ii(x^(so-to + t),t) for tG(-oo,to], 

(39) 



Now we give a brief derivation of the integral formula for the displacement 
area discussed in Section 3.1.2. Writing x'(t; sq, to) = £^(t; Sq, to), r'(t; Iq, to) = 
ef(t;/o,to) and a'(t;/o,to) = £^a(t; /q, to), (15b) takes the form: 

a(t; /o, to) = r(t; Iq, to) |u(x'^(s)) | = u(x'^(so - to + t)) A x(t; sq, to) , 

(40) 

Differentiating (40) with respect to r gives 

-^a{T; Iq, to) = trace{Z}xn(x'^(s - t + t))} a{T; Iq, to) + fJ,'^{s -t + T,T) , 
dT 

(41) 



that has a clear resemblance to (10). Unlike (10) that is two-dimensional for 
x'(t; So? ^o), (41) is a scalar, linear ordinary differential equation and hence can 
be solved analytically. Doing so with initial condition a'(to;lo,to) = gives 
(16). 



A. 2 Velocity Field Defined by a Data Set 

When the velocity field u(x, t) is given as a data set, we need to consider more 
carefully the applicability of the collection of mathematical results discussed 
above since the data set does not contain an explicit e. Practically, this means 
we must first decide on a particular decomposition of u(x, t) into il(x, t) and 
eu(x, t). Although the choice of u(x) is not unique, a natural choice is the 
time average (mean) of u(x, t). The geometrical structure of the fiow defined 
by u(x), we can make a choice for C = {x'^(s)}. The residual becomes the 
unsteady fiuctuation, i.e., 

u'(x, t) = u(x, t) - u(x, t) . (42) 

Once a decomposition of the data set is chosen we then need to examine the 
smallness of u'(x, t) with respect to u(x). As in Section 3 and Appendix A.l, 
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we do so for the separate types of C: finite time interval, and semi- or bi-infinite 
time interval. 

The essential result required is to show that the unsteady perturbation is 
actually small for a finite-time interval C. There are numerous ways of doing 
this. One way is to consider the quantity: 



where | ■ | denotes a convenient norm. This quantity can be estimated numer- 
ically in a neighborhood of C for a time interval of interest. 

For either a semi-infinite or a bi-infinite time interval, a'" (x, t) defined by (43) 
will not generally be small at the endpoints of C where we have x'^''* since 
|il(x'^'^*)| = 0. However, this does not affect the validity of the TIME functions 
since the integrand of each function contains il(x'^(so — ^o + ^)) which vanishes 
exponentially fast as x'^*^* is approached along C . The argument here is the 
same as the proof of absolute convergence of the Melnikov integrals given in 
the references at the beginning of this appendix. Effectively, the smallness of 
(43) only needs to be established along C outside of a neighborhood of any 
endpoints of C that are hyperbolic stagnation points. 



B Relation to Lagrangian Transport 



This appendix describes the similarities and differences between the Lagrangian 
lobe dynamics that is closely associated with the pseudo-lobe dynamics of the 
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TIME method along W . All TIME functions used in this section have the 
superscript {•}'^, reflecting the fact that the Eulerian boundary C we are 

TT 

considering is W . We begin by presenting a brief description of Lagrangian 
transport from the dynamical systems point of view; details of the theory can 
be found in [11,23,18,2]. 



B. 1 Overview of Lagrangian transport 



Geometry of Lagrangian invariant manifolds. Lagrangian lobes are formed by 
time-dependent unstable and stable invariant manifolds, VV'^(t) and VV'^(t), 
of upstream and downstream time-dependent DHTs, xl'^(t) and x^*(t). The 
geometry is schematically shown in Figure 8. For convenience, we parameterize 
the coordinate of the points on W^it) by x^(s,t) and on >V'^(t) by x'^(s,t). 
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respectively, by the flight-time coordinate s using their normal projection onto 
W" = {x^(s)}, i.e., (x^'^(s,t) - x^(s)) ■ ^f (/^(s)) = for both W^(t) and 
yV^{t). Then the normal distance from W to W'^(t) and VV'^(t) is: 



^U,Sf,Hf^ +^_(^U,S(, +\_^Hf^^ . ^CfiC 



rrr(5),t) = (x-'^(s,t)-x-(.))A^^(^(s)) , (44) 

where the subscript {■}± here represents the normal projection. Accordingly, 
the normal distance from VV'^(t) to W^(t) is defined by: 



ri(/,t)=r^(/,t)-ri(/,t). (45a) 

The superscript {■}^ stands for Lagrangian. This normal distance r^{l,t) is 
approximated using the so-called Melnikov function up to leading order, which 
is identical to the displacement distance function of TIME: 



r^r(«),t) = SM = r^H^'islt). (45b) 



Principal Intersection Point (PIP) sequence. An intersection sequence between 
yV^{t) and W'^(t) can be expressed as a discrete sequence, {Sj(t)} with 
sfit) < Sj^^(t). These intersecting points are called principal intersection 
points (PIPs). The sequence can be transformed to the arc-length coordi- 
nates {Ifit)} using lf{t) = l^{s^{t)). Invariance of W^(t) and W^{t) guar- 
antees that a trajectory starting at any PIP will remain a PIP for all time. 
From (45b), the PIP sequence can be approximated by the zero sequence of 
r^{l^{s),t) up to leading order: 



sf{t) - sf{t), lf{t) ~ P(sf (t)). (46) 



Lagrangian lobe and its classification by the directionality of transport. A La- 
grangian lobe, {Cjj_^^^{t)}, is then defined by segments of W^(t) and W^(t) 
between a pair of adjacent PIPs corresponding to Sj'it) and Sj'_^_i(t). We classify 
them into the two types, C^^ or C^'^ , depending on whether the corresponding 
segment of W'^(t) lies to the left or to the right of the corresponding segment 
of W'^(t), where the directionality is measured with respect to the forward 
direction of s along W . 

Area of a lobe. Using an elaborate computational scheme, it is possible to 
compute the area A{Cjj_^_-^^(t)) of Cjj_^_i{t) [15]. Using (45b) and (46) as well 
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[Fig.8] 



as (29) on W , leading order approximation of A{Cj.j_^_l(t)) is given by the 
Melnikov function [11], which we denote by A(£j"-^^(t)): 



^(4.+iW) ^ ^(4.+iW) =^(^£+iW) ' (47) 

where A{Lfj^^{t)) is given by (29). 

Moving boundary and boundary PIP. In order to describe Lagrangian trans- 
port, a Lagrangian boundary must be defined using the segments of moving 
invariant manifolds VV^(t) and W'^(t): 

W\t) = {x^(s, t), s < s\t)} U {x^(s, t), s > s\t)}, (48a) 

where s^{t) C {Sj{t)} is called a boundary PIP as schematically shown in 
Figure 8a. Unlike Eulerian transport with a well-defined stationary boundary 
W , the selection of s^{t) out of all the existing PIPs and hence the selection 
of VV^(t) may not be unique. However a physically meaningful choice may be 
to have VV^(t) = {x^(s,t)} geometrically "close" to W [10]. This condition 
requires smallness of both \r^{l^{s),t)\ and \r']_{l^{s),t)\ in the neighborhood 
of s^(t), and hence smallness of |^i(/'^(s), t)| approximated by |r'^(/^(s),t)|. 
From (45b), a reasonable choice of s''{t) is therefore 



s\t) = {se s^{t) I min\l"{su) - l"{s)\} (48b) 



H 



where Su = {s | 7T2aa;|u(x'^(s))|} is the maximum velocity point on W 

Further classification of C^'j^i^t) by the timing of transport. The relation be- 
tween s^(t) and Sj'{t) leads to a further classification of the Lagrangian lobes 
beyond C^^ or C that also describes the timing of transport associated with 



Particles in Cjj_^i{t) with Sj{t) < s^{t) are yet to cross VV^(t) and hence 
have not involved in the transport process yet. Because where VV^(t) is a 
segment of W^(t), we classify these Lagrangian lobes as C(^=)'^^ and C^^(='') 
with the same directionality as C^^ or £'^^, respectively On the contrary, 
particles in Cj'j_^_^(t) with Sj(t) > s^(t) have crossed W^(t) and hence have 
already involved in the transport process. We classify these Lagrangian lobes as 
£US{=b) g^j^^ £{b=)su ^\^\^ ^]-^g same directionality as C^^ or C^^ ^ respectively. 

Discrete transport time sequence. As particles advect downstream in the un- 
steady flow, s^{t) must change in time not only by advection but also for 
keeping the shape of W^{t) as close as possible to W as defined in (48b). 



27 



Accordingly, s^(t) must jump from Sj(t) to sj'_^ (t) at some discrete time se- 
quence, {t\.} where jl is the corresponding PIP number and Ajl and increment 
for the A;-th jump. The integer Aj^ should be positive so that the shift of s^{t) 
moves towards upstream at t^. Between two consecutive jumps t G (t^,t^_,_i), 
there is no Lagrangian transport between the kinematically-distinct regions 
because they are defined by the same (moving) boundary >V^(t). At the dis- 
crete time sequence {t^.}, Lagrangian transport by the lobe dynamics occur 
due to the change of the boundary. 

Turnstile mechanism of transport. At the discrete time sequence {t^}, a total of 
Ajt-1 lobes, £^,„+i(4), for n = Ut-Ajl), ..., (jt-1), go across W\tl) into 
the other region. Therefore, the lobes of £(^=)'S'^ turn into C^^^=^) and those 
of C(^=)'^^ turn into £'5'{/(=fe)_ ^jj-^jg jg ^j-^g so-called turnstile mechanism. As a 
consequence, Lagrangian transport involves only a few lobes at a time that 
are spatially confined near the time-dependent boundary point (s^(tfc),t^). 



B.2 Pseudo-lobe dynamics and Lagrangian lobe dynamics 



The turnstile mechanism of the Lagrangian method is crucial in understand- 
ing the transport of particles between the kinematically-distinct regions across 
the invariant boundary VV^(t). Lagrangian transport occurs locally only near 
x^(s^(t^), t^) as above, at the discrete time sequence {tl}. Therefore each 
Lagrangian lobe Cjj_^_i{t) can be classified decisively whether transport asso- 
ciated with it has occurred or not, by its relative position with respect to the 
s\t). 

Table 3 presents a general comparison between the TIME and the Lagrangian 
method. The TIME method computes the net amount of transport over time 
across a stationary Eulerian boundary C. At time t, pseudo- lobes {L^^_,_j^(t)} 
may contain a mixture of transport that occurred in the past and in the future. 
This effect particularly manifests itself along the heterolinic connection W , 
for the past {—oc,t] and the future r G [t, cxo). This is because the transport 
boundary is fixed as W unlike VV*(t), which changes from VV^(t) to W'^it) at 
s'\t). Therefore the TIME method by itself does not present such a mechanism 
because it deals with the net amount of transport over time but not individual 
particles. 

A parallel development of the turnstile mechanism for the TIME method is, 
however, possible by recognizing the nature of the TIME method and realizing 
the meaning of the Lagrangian transport boundary VV^(t). It requires the 
following modification to R^{t) as follows, where R^{t) is defined in (25) along 

TT TT 

W . In the downstream direction where transport has happened and W 
corresponds to W^it) associated with the downstream DHT, W^it) should 
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correspond to R^{t). In the upstream direction where transport is yet to 
happen and W corresponds to W^{t) associated with the upstream DHT, 
however, W"^(t) should correspond to the mirror image of R^{t): 



f?®(t) = {(/,r)|r = -r^(/,t)} (49) 

see Figure 9, in comparison with Figure 2. Table 4 gives the summary of the 
definitions below. [Fig. 9] 



Pseudo-manifolds. We define the unstable pseudo-manifold W'^{t) and stable 
pseudo-manifold W^(t) by the segments of W , R^ (t), and -R®(t) as follows. 
In the arc-length coordinate, W^(t) = {(/,r) | (/,r^(/,t))} and W^{t) = 
{(/,r) I (/,r'^(/, t))} are defined by 



0, for / < /^(s®(t)) on W 

[ r^(/,t), for / >" (s®(t)) on R"{t) 

f -r^(l,t), for / < P(s®(t)) on i?®(t) 
[O, for/>/^(s^(t)) only 

where definition of s®(t) is given below. The distance r^(t) — r^{t) between 
py^(t) and W^{t) is r^{t) by these definitions. It also corresponds to r^{l^{s),t) 
of the Lagrangian transport (45b). 

Pseudo-boundary PIP. We select s^{t) to be 



S^{t)r^s\t) (51) 

out of all the existing {s^(t)} so that the resulting ly^(t) and W^(t) are 
geometrically closest to their corresponding Lagrangian invariant manifolds 
W^(t) and W^it). This choice of s®(t) implies large |u(x^(s^(t)))| from (48b) 
and hence small |r^(/^(s®(t)), t)|. 



Pseudo-Lagrangian lobes. Originally the pseudo-lobes are defined for C in Sec- 
tion 4.2 by the segments of R^{t) and C. For transport dynamics across W , 
however, we follow the convention of the Lagrangian lobes defined in Sec- 
tion B.l. The pseudo-Lagrangian lobe is defined by segments of W'^{t) and 
W^{t) between a pair of adjacent pseudo-PIPs corresponding to s^{t) and 
sf^iit). We classify them into two types, L^*^ or L"^^, depending on whether 
the corresponding segment of W^{t) lies to the left or to the right of the cor- 
responding segment of W^{t) in the forward direction of s along the boundary 

TT 

W . Table 2 was computed by following this convention. 
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[Tab. 4] 



Directionality of transport across W by the pseudo-Lagrangian lobes. Like La- 
grangian lobe dynamics, the relative position of s^{t) with respect to s®{t) 
determines whether or not the j-th pseudo-Lagrangian lobe has been trans- 
ported across W . This leads to a further classification of pseudo-Lagrangian 
lobes, L(^=)^^ and L^^(=^) for L^^ as well a L(^=)^^ and L^^(=^) for L^^ . 
The pseudo-Lagrangian lobes of L^^^^^^ and L^'^(=^'> lie before s®(t) and 

TT 

will cross W in the future r > t. In contrast, the pseudo-Lagrangian lobes of 
]JJS(=H) g^^j j^(H=)su jjg after s®(t) and have crossed W in the past r < t. 

Discrete transport time sequence of TIME. As in the case of the boundary PIP 
s*(t), the inverse pseudo-PIPs s'^it) also forms a discrete time sequence {tf }. 
By the construction of s®(t) (51), {tf } is an approximation of the discrete 
time sequence of the boundary PIPs {t\}, i.e., tf ~ t^ up to leading order. 

Turnstile mechanism of the pseudo-Lagrangian lobes. Like the Lagrangian lobe 
dynamics, at the discrete time sequence {t^}, a total of Ajf — 1 pseudo- 

TT 

Lagrangian lobes go across W into the other region. Therefore, the lobes of 
l(H=)su ^^^^ jj^^Q iUS(=H) ^^^ ^Yiose of 1^^=)^^ turn into 1^^^=^) . This is 

the so-called pseudo-turnstile mechanism of the TIME method. 
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Lagrangian lobes using the pseudo-manifold described in 
Appendix B and by the Lagrangian lobes using the method 
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Results of the TIME pseudo-Lagrangian lobes are based on 
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area of a pseudo-lobe and a Lagrangian lobe corresponds to 
the amount of transport: a positive value represents transport 
from the subpolar gyre to the subtropical gyre, while a 
negative value represents transport from the subtropical gyre 
to the subpolar gyre. 
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Lagrangian method. 

Relation between the TIME method and Lagrangian lobes for 

TT 

transport associated with W . For (pseudo-)lobe type, "right 
to left" and "left to right" indicate the direction of transport 
while "future" and "past" indicate the timing of transport 
with respect to the present time t. See text for details. 
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List of Figures 



Geometry associated with tlie TIME metliod: a) Eulerian 
boundary C = {x'"(s)} and tlie mean velocity 0verlineu{'x); 
b) instantaneous flux 0(x, t) across C as tlie parallelogram in 
the unsteady flow; c) accumulation mP{s,t]D) (the shaded 
portion represents the accumulation); d) displacement distance 
mP{s,t]D) (the shaded area represents a positive pseudo-lobe 
area). HOl 



Schematics of pseudo-lobes along C\ a) in x; b) in (/,r); c) in 
{s,m). The shaded areas correspond to the net accumulation 
of fluid particles originally on the right side of C, where the 
direction is deflned with respect to the direction of increasing 
s, which accumulated onto the left side of C at time t. 

Double-gyre circulation: a) streamfunction ipi^Xjt) at t = ^944; 
b) unstable manifold W'^(t) (solid line) and stable manifold 
W^it) (dashed hne) of the inter-gyre transport by the 
Lagrangian method; and c) reference streamfunction ip^x) 

and the Eulerian boundary W = x'^(s) of the inter-gyre 
transport by the TIME method. In (a) and (c), contour 
interval 2000 with dash lines for negative values. In (b), 
numbers correspond to those in Table 2 In (c), PF = x^(s) is 
the thick solid line with four white circles plotted at every 250 



days starting from in s 


= 0, 


i.e. 


, x^(^S) 


5 X (S250), X (S500), 


and x'^(s^Qoo); the diamonds at 


x^(sj), 


x^(sn), and x^(ss) 


with (sj,sn,ss) = (s^io, 


"^129 


) "^174.5)' x ( 


Sj) almost overlaps 


x^(.S). 











The instantaneous flux function fi^{s,t) in the HovmoUer 
diagram with dash-dot line for the negative values and 
contour interval 40km^/day. The diagonal line is an example 
of reference trajectory going through (so,to) = (0,1300). 
The horizontal lines corresponds to the four phases during 
one period of the oceanic oscillation, at t = ^944, ^944 + T/4, 
t;44 + T/2,t;44 + 3T/4. 

The instantaneous flux function /i^(s, t) at t = tg44, ^944 + T/4, 
^944 + T/2, t;44 + 3r/4 from bottom to top (right) with 
corresponding time indicated in Figure 4. 
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6 The accumulation functions a^{s, t) (dash hne), a^{s, t; t) (sohd 
hne), and a^{s,t;t) (dash-dot hne) at ^944, tgg2(= ^944 + T/4), 
^*io2o(= ^944 + T/2), and t*058(= ^944 + 3T/4) (from bottom to 
top) with corresponding time indicated in Figure 4: a^{s,t) 
almost overlaps a^{s, t; t) for s < sj and a^{s, t; t) for s > ss- US] 

7 The displacement distance functions, r^{l,t) (dash line), 
r^{l,t]t) (solid line), and r^{l,t;t) (dash-dot line) at ^944 
corresponding to the bottom panel in Figure 6; r^{l,t) almost 
overlaps r^{l,t]t) for / < l^{sj) and r^{s,t;t) for / < /^(ss), 
although pseudo-lobes of r-'^(/,t) and r^{l,t]t) accumulate 

near / = and cannot be seen in this figure. HE] 

8 Schematics of Lagrangian lobes in x. The shaded lobe 
Cj_^2~j+?,{^) corresponds to fluid particles which have already 
transported from the right to the left of the deformable 
Lagrangian boundary >V^(t), while the swatched lobe 
^jj+i {t) corresponds to fluid particles to be transported in 

the future time. HT] 

9 Schematics of pseudo-Lagrangian lobe in a) x and b) (/, r). The 
shaded pseudo-Lagrangian lobe Lj^2~j+3i't) corresponds to net 
amount of fluid particles which have already been transported 
from the right to the left of W , while the swatched lobe 
Ljj+i (t) corresponds to the net amount of fluid particles to 
be transported in the future time. 
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Symbol 


Definition 


Equations 


s 


flight-time coordinate along C 


(4) 


C = {x^(s)} 


Eulerian boundary over the segment s £ [sa,si,], including: 


(4) 




W semi-infinite for s G [so,oo) 


(6) 




W semi-infinite for s G ( — 00, sq] 


(6) 




TT 

W bi-infinite for s £ (—00, 00) 


(7) 


{s,t) = {so-to + t,t) 


reference particle trajectory (advection) along C 


(5) 


{l,r) = {l^\s),r) 


arc-length coordinate in two dimensions 


(8) 


wP{s,t;tQ : ti) 


accumulation function 


(14),(19)-(21) 


a^{s,t;to:ti) 


displacement area function 


(15),(19)-(21) 


r^'{s,t;to:ti) 


displacement distance function 


(15),(19)-(21) 


sf(t) 


pseudo-PIP 


(26) 


Lf.^iit) 


pseudo-lobe, which belongs to the following two types: 
]^RC fpom right to left across C 
j^CR fj-om left to right across C 


(28) 


^(^^.+iW) 


signed area of pseudo-lobe 


(29) 



Table 1 

Glossary of TIME functions. See text for details. 
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TIME 


Lagrangian 


lobe no 
id. ij) 


(km^) 


(km2) 


1 


6858 


6258 


2 


-6776 


-6355 


3 


6796 


6201 


4 


-6783 


-6216 


5 


6760 


6424 


6 


-6838 


-6515 


7 


6819 


6696 


8 


-6811 


-6119 



Table 2 

Quantitative comparison of transport by the TIME pseudo-Lagrangian lobes using 
the pseudo-manifold described in Appendix B and by the Lagrangian lobes using 
the method described in [15] at t^^^ (see Figure 3b for the lobe number). Results of 
the TIME pseudo-Lagrangian lobes are based on the pseudo-manifolds described in 
Appendix B. The signed area of a pseudo-lobe and a Lagrangian lobe corresponds 
to the amount of transport: a positive value represents transport from the subpolar 
gyre to the subtropical gyre, while a negative value represents transport from the 
subtropical gyre to the subpolar gyre. 
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TIME 


Lagrangian 


type 


particle, flow property Q 


particle 


boundary 


spatial 


any reasonable C 


W\t) only 


temporal 


stationary 


continuously deformable 
& discontinuous at {t^} 


u(x,t) 


restriction 


u' (x, t) small along C 


no restriction 


spatial requirement 


along C only 


non-local in x 


temporal period 


flexible 


inflnite 


transport 
process 


temporal 


continues 


turnstile at {t^} 


spatial 


all along C 


locally near x''(i* (4), t^) 


computation 


efficient 


exhaustive 


accuracy 


leading order 


accurate 



Table 3 

A general comparison between the TIME method and Lagrangian method. 
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TIME 


Lagrangian 




intersection sequence 


{sfm 


{sfm 




geometry 


distance 


r^{l,t) 


r^il,t) 






area 


^(^f.+i(*)) 


A{Cf^^^,{t)) ^ A{Cf^^ 


+1W) 




boundary point 


s^it) 


s\t) 




transport by 


time sequence 


{tf} 


{4} 




turnstile 


lobe number 


ft 


Jl 






lobe increment 


Ajf 


Ait 




(pseudo-) 


right 


future 


L(h=)us 


^{b=)US 




lobe 
type 


to left 


past 


lUS{=H) 


^US{=b) 




left 


future 


lSU(=H) 


^SU(=b) 






to right 


past 


L(h=)su 


^{b=)SU 





Table 4 

Relation between the TIME method and Lagrangian lobes for transport associated 
with W . For (pseudo-)lobe type, "right to left" and "left to right" indicate the 
direction of transport while "future" and "past" indicate the timing of transport 
with respect to the present time t. See text for details. 
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b) u'(x,t) 





c 

rc{lc(s),t;D) 



-4- ^ 



Fig. 1. Geometry associated with the TIME method: a) Eulerian boundary 
C = {x*^(s)} and the mean velocity 0fer/ineu(x); b) instantaneous flux (/>(x, t) 
across C as the parallelogram in the unsteady flow; c) accumulation m {s, t; D) (the 
shaded portion represents the accumulation); d) displacement distance mP{s,t]D) 
(the shaded area represents a positive pseudo-lobe area). 
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a) 







l',Jt} 



b) 





,W 



Fig. 2. Schematics of pseudo-lobes along C: a) in x; b) in (/,r); c) in (s,m). The 
shaded areas correspond to the net accumulation of fluid particles originally on 
the right side of C, where the direction is defined with respect to the direction of 
increasing s, which accumulated onto the left side of C at time t. 
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a) ip(x,t) 



b) w"(t) & W^(t) 



l/H : 



c) W in Reference Flow 




I, 1000 





200 400 600 800 1000 
X (km) 



200 400 600 800 1000 200 400 600 800 1000 

X (km) X (km) 



Fig. 3. Double-gyre circulation: a) streamfunction ■ip(x,t) at t = t^^^; b) unstable 
manifold W^(t) (solid line) and stable manifold W*(t) (dashed line) of the inter-gyre 
transport by the Lagrangian method; and c) reference streamfunction iIj{x) and the 
Eulerian boundary W = x (s) of the inter-gyre transport by the TIME method. 
In (a) and (c), contour interval 2000 with dash lines for negative values. In (b), 
numbers correspond to those in Table 2 In (c), M^ = x^(s) is the thick solid line 
with four white circles plotted at every 250 days starting from in s = 0, i.e., x (sq), 
x^(s25o), ^^('*50o)' ^^'^ ^^('^looo)' ^^^ diamonds at x^(sj), x^(sn), and x.^{ss), 



with (sj,SN,ss) = (•5*10' ^129) •5*74.5); ^ (sj) almost overlaps x^ {sq) 
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1500- 



1400 



1300 



1200 



1100 



1000- 



900 




500 



Fig. 4. The instantaneous flux function fj, {s,t) in the HovmoUer diagram with 
dash-dot line for the negative values and contour interval 40km^/day. The diagonal 
line is an example of reference trajectory going through (so,to) = (0,1300). The 
horizontal lines corresponds to the four phases during one period of the oceanic 
oscillation, at t = ^944, t944 + T/4, ^944 + T/2, tl^^ + ST/A. 
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Fig. 5. The instantaneous flux function /i^(s,t) at t = tg^4^, ^944 + T/4, tg44 + T/2, 
^944+3r/4 from bottom to top (right) with corresponding time indicated in Figure 4. 
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Fig. 6. The accumulation functions a {s,t) (dash hne), a {s,t]t) (sohd hne), and 
a^{s,t;t) (dash-dot hne) at ^§44, t982(= ^944 + T/4), tio2o(= *944 + T/2), and 
*1058(~ ^944 + 3T/4) (from bottom to top) with corresponding time indicated in 
Figure 4: a^{s,t) almost overlaps a^{s,t;t) for s < sj and a^{s,t;t) for s > ss- 
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Fig. 7. The displacement distance functions, r {l,t) (dash hne), r {l,t;t) (sohd 
Une), and r^{l, t; t) (dash-dot hne) at tg^^ corresponding to the bottom panel in Fig- 
ure 6; r^{l,t) almost overlaps r^{l,t;t) for / < l^{sj) and r^{s,t;t) for I < l^{ss), 



although pseudo-lobes of r^{l,t) 
be seen in this figure. 



and r^ {l,t;t) accumulate near / = and cannot 
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Fig. 8. Schematics of Lagrangian lobes in x. The shaded lobe >C -.g • , 3(i) corre- 
sponds to fluid particles which have already transported from the right to the left 
of the deformable Lagrangian boundary W^(t), while the swatched lobe L- Z_\ (t) 
corresponds to fluid particles to be transported in the future time. 
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Fig. 9. Schematics of pseudo-Lagrangian lobe in a) x and b) {l-,r). The shaded 
pseudo-Lagrangian lobe L-_^2-_^j^{t) corresponds to net amount of fluid particles 



H 



which have already been transported from the right to the left of W , while the 



swatched lobe L- 



{H=)US 



ij+1 



(t) corresponds to the net amount of fluid particles to be 



transported in the future time. 
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